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1. Introduction 

In equilibrium statistical mechanics, one can prove that the grand canonical free energy of the 
system can be written as a functional of the one-body density only f\\. The density distribution 
which minimizes the grand potential functional is the equilibrium density distribution. This 
statement is the basis of the equilibrium density functional theory for classical fluids which 
has been used with great success to describe a variety of inhomogeneous fluid phenomena 
Q. Out of equilibrium there is no such rigorous principle. However, macroscopically one 
can find phenomenological equations for the time evolution which are based on macroscopic 
quantities only, e.g. the diffusion equation, the heat transport equation and the Navier-Stokes 
equations for hydrodynamics. 

In perhaps one of the most simple microscopic cases, a system of interacting Brownian 
particles, a number of equations for the time evolution of the particle density have been 
proposed which aim to incorporate results from equilibrium density functional theory. These 
equations are usually referred to as dynamical density functional theory (DDFT) or time 
dependent density functional theory. Most of the equations take a rather similar form but 
there has been some controversy as to whether the time evolution ought to be stochastic 
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10 in 13 1^1 171 or deterministic (HI HOI in nature. In reference Q the authors attempt to 
address these issues. However, we believe further clarifications are necessary. The reason 
for the controversy is that one has to state precisely what type of density distribution one is 
talking about, namely whether one speaks of the instantaneous density p{r,t) at time t, an 
ensemble averaged density p{r,t), or a spatially and/or temporally coarse grained density 
p{r, t). We will give proper definitions for these densities later. The time evolution equation 
for p{r, t) is discussed in references lEHBElISll and for Brownian particles this is, of course, a 
stochastic equation. p(r, t) results from taking an ensemble average over the stochastic noise 
and is therefore uniquely valued at time t; the equation governing its dynamics must therefore 
be deterministic - this is the approach in references |l8l|9l[T0||. The equation governing the 
dynamics of p(r, t), the coarse grained density, will, of course, still contain a noise term. The 
main objective of this paper is to clarify the relationships between these three approaches and 
to discuss the similarities and differences. In order to visualize better the connections between 
the DDFTs for /)(r, t), p(r, t) and p(r, t), we employ the flow chart in figure[l] 

This paper proceeds as follows: In section |2l we formulate the dynamics of a system of 
interacting Brownian particles. We derive a deterministic DDFT for the ensemble (or noise) 
averaged density in section |3l In section |4l we formulate the stochastic DDFT for the coarse 
grained density and finally, in section |5l we discuss the results. 

2. Brownian dynamics of interacting particles 

The starting point for our considerations is the Langevin equation for a system of N identical 
Brownian particles with positions {T'i(t)}j=i...Ar (see the dashed box in figure[T]). The particles 
interact pairwise via the (pair) potential V{r). In addition, there is a one-body external 
potential $(rj, t). This equation reads: 



where F{ri — Vj) = —'ViV{\ri — rj\) is the interaction force between a pair of particles and 
G{ri, t) = — Vj$(rj, t) is the external force on particle i. The uncorrected Gaussian random 
force on each particle {T7i(t)}j=i...Ar is fully characterized by 



The particles are labelled with lower roman indices and the spatial directions (i.e., x, y, and z 
in three dimensions) are labelled with upper indices. A possible prefactor (mobility constant) 
in front of the deterministic part of ([T]) can be eliminated via the Einstein relation and a 
corresponding rescaling of time. T is the temperature. Note that the noise term in dU) is 
additive, i.e. it does not have a prefactor which depends on the particle position. In this case, 
the Ito and the Stratonovich stochastic calculus are equivalent ifTDl . 

In the next step, we define the instantaneous microscopic density p{r,t;{ri(t)}) = 
J2iLi — f)- This density depends implicitly on the actual positions of all particles 

{rj(t)} and it should be considered a density 'operator' rather than a measurable particle 




(1) 



j=l...N 



{mit))=0 and (r/f(t)?7-(t')) =2T5,,y™5(t-0. 



(2) 
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density. Using Ito stochastic calculus, one can show that the time evolution of this sum of 
delta functions is given by (see references lEl equation (8)] or |l5l equation (17)]): 



dt 



p{r, t) G{r, t) - p{r, t) J p{r' , t) F{r - r') dV 



-TV/)(r,t) + v//3(r,t)^(r,t) 



(3) 



i.e., an equation which, at a first glance, does not depend upon the individual particle positions. 
This step corresponds to path 2, which leads to the left dotted box in figure [T] The random 
forces on the particles rii{t) have been replaced by a Gaussian random field ^{r,t) with 

(|(r,t)) = and {^^{r,t) ^""{r' ,t')) = 2T 6^"" d{r - r') 6{t - t'). (4) 

Introducing the 'energy' functional 7Y[/)] iQ 



n[p] = j |Tp(r,t) [lnA3p(r,t)-l] +p(r,t)<l>(r,t) 
+ \ j p{r^)V{\r-r'\)p{r\t) dVj dV, 



where A is the de Broglie wavelength, Q can be written as 



dp{r,t) 
dt 



p{r,t)V 



(5) 



(6) 



6p{r,t) 

In contrast to the noise in equation ([T]), the noise in Q is multiplicative. However, as shown 
in ifTll . for a conserved multiplicative noise, Ito and Stratonovich calculus are equivalent and 
therefore Q, the Langevin equation for p{r, t), should have the same form when obtained by 
Stratonovich calculus. 

We denote the probability density for finding the particles in the system at time t with 
the positions {re}i=i,„N by w{ri, . . . , tm, t). The time evolution of this probability is given 
by the Fokker-Planck equation for dU), namely 

dw{{ri},t) 



dt 



G(r„t)-TV, 



w{{re},t). 



(7) 



1=1. ..iV \_j=l...N 

In figure [T] the derivation of ^ is denoted by path 1 . Equation Q is also referred to as 
the Smoluchowski equation, namely the drift-diffusion equation for a many particle system 

cniioi. 

However, we are more interested in the particle density distribution, rather than in the 
positions of the individual particles. The probability of finding the system in a configuration 
with the density p{r, t) is given by 



W[p] 



wi{n}, t) S [p(r, t) - p{r, t- {r,})] dVi . . . dV^v- 



(8) 



p{r,t; {vi}) is a sum of (5-functions and therefore, from the above definition, W[p] can 
only be non-zero when p(r, t) is also a sum of exactly A^ ^-functions. In order to make a clear 
distinction between smooth densities and spiky sums of delta functions we will keep a hat on 
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the argument of W. The time evolution for W[p], which in fact is the Fokker-Planck equation 
corresponding to is given by fJ] (see also iflBl or appendix 8 in lfT4l ): 



dW[p] _ 
dt ~ ~ J 6p{r,t) 



V-/5(r,t) V 



■y. S , 6n[p] 



W[p] \ dV (9) 



5p{r,t) 5p{r,t)_ 

with the functional Ti.[p\ as given in (jSj). The derivation of Q from the Langevin equation 
® is path 5 connecting the two dotted boxes in figure [T] Equivalently, one could also obtain 
^ via a change of variables from {r-jjj^i jv in ^ to p(r, t); this is path 9 in figure[T] The 
equilibrium distribution corresponding to ^ is given by VF'^^'')[p] oc exp [—^ ('^[p] + p p)] 
with an arbitrary chemical potential p. The freedom of choice for p comes from the fact that 
^ has the form of a conservation law. 

Eq. ^ together with equation Q is a DDFT for p{r, t), i.e. an equation of motion for 
the instantaneous microscopic density. However, although equation ^ is written in terms 
of p{r,t), this sum of delta functions contains the exact positions of all particles (modulo 
permutations) and the solution of © is formally equivalent to integrating O - see also the 
appendix in reference ifTSl . Generally, one is not interested in the individual trajectories 
of the particles, rather in the time evolution of quantities (in particular the fluid density) 
averaged over many realizations of the noise or, alternatively, in the evolution of coarse 
grained quantities. 



3. Ensemble averaged density 

We define the ensemble averaged density as the density p{r, t) averaged over many solutions 
of Q with equal initial conditions but different realizations of the noise term ^(r, t), i.e. 

p(r, t) = J pir, t; {r,}) w{{r,}, t) dVi . . . = (p(r, t)). (10) 

In other words, in order to measure p(r, t) one solves © many times for different 
realizations of the noise ^(r, t) and averages the result at each time t. Since the particles 
are indistinguishable, w{{ri},t) is symmetric with respect to permutations of particles and 
therefore p(r, t) = N J w{r, r2, . . . , tn-, t) d?r2 . . . d?rN. We can obtain the time evolution 
of p(r, t) by integrating out — 1 degrees of freedom in Q to obtain lfT0lfT6ll : 

= - V ■ I [G(r, t)-TV] p(r, t) + j p^^) (r, r', t) F{r - r') dvj , (1 1) 

where 

p(2)(r, r', t) = N{N-l) j w{r, r\ rg, . . . , r^v, t) dVg . . . dV^, (12) 

is the two-body density distribution function. We have assumed that surface terms which 
appear in the partial integrations vanish. Here we see that for interacting particles, the time 
evolution of the one-body density distribution p{r,t) depends upon the two-body density 
distribution function, i.e. on the particle correlations in the system. In fact, (fTTT) is the first 
equation of a hierarchy of equations for the n-body density distribution functions (similar to 
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the BBGKY hierarchy) which one obtains by integrating out (A^ — n) degrees of freedom in 
© -see Col- 

So far, everything is exact, but in order to make use of this hierarchy one needs a closure. 
One way, which is employed in reference lITOll . makes use of a result from equilibrium density 
functional theory, in which there is an exact sum rule connecting the equilibrium two-body 

('2) 

density distribution function peq{r,r') to the gradient of the one-body direct correlation 
function c^^^ (r) 

kBTp{r) Vc^^\r) = J pfj{r,r') F{r - r') dV. (13) 

In equilibrium c^^^(r) is the effective one-body potential due to interactions in the fluid and is 
given by the functional derivative of Jix [p] > the excess (over ideal) part of the Helmholtz free 
energy functional, —kBTc^^\r) = IDIlQl, and we obtain the following exact result 

p{r) = - j pg) (r, r') F{r - r') dV. (14) 

Assuming this relation remains valid for the non-equilibrium fluid is equivalent to assuming 
that the two-body correlations in the non-equilibrium fluid are the same as those in an 
equilibrium fluid with the same one body density profile |l8l[T0|. Within this approximation 
we can write the time evolution equation (fTTT) as lISlfTOl: 

dp{r,t) 



dt 



(15) 



5p{r,t) 

with the Helmholtz free energy functional! 

rW,t)] = j p{r,t) [T (lnAV(r',t)-l)+$(r,t)] dV + ^,[p(r, t)]. (16) 

Since p{r, t) is a quantity averaged over the realizations of the noise, the time evolution 
equation ( fT5l is necessarily deterministic. Eq. ( fT5l was originally written down (without 
derivation) by Evans [1|. The derivation of dTSt presented in this section corresponds to 
path 3 in figure[T] In reference JHl Marconi and Tarazona follow path 4 in figure[l]and ([131) is 
derived from Q by performing an ensemble average. Note that dTSt also remains valid when 
the particles interact via multi-body potentials ifTDIl . 



4. Coarse grained density 

The ensemble averaged density p{r, t) discussed in section|3]is the quantity often considered 
by theorists, but it is not the quantity one would measure by performing one individual 
experiment (numerical or real). In order to measure the density as a function of time in 
an experiment one must define a probe volume in space and a time window and count the 
number of particles in that probe volume averaged over the time window. Since we intend to 
resolve the atomistic structure of the fluid, the spatial resolution has to be very good. In other 

X Note that the equilibrium excess Helmholtz free energy functional Jix [p], is only known exactly for the case 
of hard-rods in one dimension HIT] (see also 0). In practice one can only obtain an approximation to the exact 
functional. Some of the more recently derived functionals are rather sophisticated and very accurate Q. 
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words, we assume an experimental resolution function of the form Ki{r, t) = K{t) 5{r) and 
then the measured density is 

p(r,t) = j K{t-t')p{r,t')dt'. (17) 

For simplicity we assume further that K(t) has a compact support. Similarly, we can define a 
two-particle probe, with a resolution function of the form K2{r, r', t) = K{t) d{r) S{r'). The 
two-particle distribution function that we measure is (c.f. equation (A6) in |[T]|) 

J K{t - t') p{r, t') p{r', t') dt' = j Kit - t') 5{r,{t') - r)5{r,{t') - r') dt' 

+ I K{t- t') J2 ^init') - r)S{r,{t') - r') dt' 

= p^^\r, r', t) + p{r, t) 5{r' - r). (18) 

This equation defines the coarse grained two-particle density distribution function, 
pi'^)(^r, r', t). When the fluid is in equilibrium, assuming ergodicity, and if the time window 
of K(t) is large as compared to the time the system needs to explore phase space locally, 
then p^^^(r, r') ^ p^'^^r, r')\ we therefore assume p^'^\r, r', t) p^'^\r, r', t). Here, locally 
means on the scale of the correlation length and therefore this may not be true for long ranged 
correlations. 

Using Q together with (flTb and dTSt . the time evolution of p{r, t) is given by 

dp{r, t) 



dt 



T V^p{r, t) - V ■ t) G{r, t) - V • ^ p^^^ (r, r', t) F{r - r') dV 

j K{t - t') ^/p{r,t')^{r, t') dt' . (19) 



The delta function term in dTSt leads to a term containing -F'(O) = —'W{r)\r=o, namely 
the force between particles when the centers are at zero separation. We assume that this is 
zero. This will certainly be the case for particles interacting via pair potentials of the form 
V{r) = Vhsir) + Vtaiiir), where Vhsir) is the hard-sphere pair potential. Equation ( [T9b is 
exact when G does not change over time or when it only changes on time scales much longer 
than the support of K{t). However, in this form equation ( [T9b is of little use, since the right 
hand side involves p(r, t) explicitly. In the following we give a number of 'hand-waving' 
arguments which allow us to approximate the right hand side in terms of p(r, t) only. These 
arguments are far from being rigorous and are by no means a proper derivation of a DDFT for 
p(r, t). The reason for giving these arguments is that they lead to a stochastic time evolution 
equation for p(r, t) which has the same form as the time evolution equation of p(r, t), i.e. (fTSt . 
but with an additional conserved noise term. In reference H Kawasaki uses a more formal 
spatial coarse graining procedure and in reference ifTSi Perez-Madrid et. al. use a mesoscopic 
approach, to obtain a similar result. 

The time integration in the last term of equation (fT9l) in principle leads to a non- 
Markovian time evolution for p{r,t). However, on time scales large as compared to the 
support of K(t), the dynamics of p{r,t) should be approximately Markovian and we write 
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the last term of equation (fT91) in the same form as in @, but with p(r, t) replaced by p(r, t), so 
that we obtain a time evolution equation identical to (fTTT) but with an additional multiplicative 
noise term: 

dp{r,t) 



V ■ I t) - T V] p(r, t) p^^) (^.^ ^.'^ t) F{r - r') dV 
-V ■^p{^)i{r,t)\. (20) 



at 



Using the local equilibrium assumption applied in section|3l we assume equation (fT4l) with p 
replaced by p and p^"^^ replaced by p^^^ applies for the non-equilibrium coarse grained density 
distributions, and we substitute the term involving p^^^ (r , r', t) in (EUt by a term involving the 
functional derivative of the excess part of the Helmholtz free energy functional, to obtain a 
stochastic version of dTSt : 

dp(r,t) 



a, - ...... „ „ , . (21) 

Note the similarity with the time evolution equation Q for the instantaneous microscopic 
density p{r,t). However, a significant difference is that it is the Helmholtz free energy 
functional, J-'[p], given by equation (fT6b with p replaced by p, that enters into equation (ISTT) ^. 
whereas 7Y[p] entering equation Q is simply the functional as defined in Unfortunately, 
in reference Q, this distinction is not made clear. In developing (ETT) we followed path 6 
in figure [T] We can obtain the Fokker-Planck equation corresponding to equation (1211) by 
discretizing it in space. Taking the continuum limit of the Fokker-Planck equation for the 
discretized Langevin equation |l3l[T2|, we obtain an equation very similar to Q: 

This calculation corresponds to path 8 in figure[ll The equilibrium distribution corresponding 
to (l22b is given by W^^"^^ [p] oc exp [—^ {^[p] + y^t p)] with an arbitrary chemical potential p. 
This result justifies a posteriori our, initially unjustified, approximation for the noise term in 
( fT^ . which we used to obtain (ETT) . 

In this section we considered the temporally coarse grained density p(r, t) because we 
wished to make a comparison with the DDFT obtained for the ensemble averaged density 
p(r, t) in the previous section, which is used to describe the time evolution of the structure of 
the system at microscopic length scales. The derivation was by no means rigorous. However, 
a spatial coarse graining of ^ can also lead to an equation of the same form as (l22b but with 
a particular (approximate) form for the Helmholtz free energy functional - see reference Q. 
This last route follows path 7 in figure [T] 

§ More strictly, we expect the free energy functional in J2H to be a coarse-grained free energy functional, ^[p]. 
We are making the approximation !F[p] ~ ^[p]- See 1191 for a recent discussion of the difference between these 
free energy functionals. 
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5. Discussion 

In section |2l and section |4| we derived stochastic DDFTs (Langevin equations) for the 
instantaneous microscopic density p(r, t), equation Q, and the coarse grained density p(r, t), 
equation d^TTl . respectively, as well as the corresponding Fokker-Planck equations. In section|3l 
we also derived a deterministic DDFT, equation dTSt . for the ensemble (or noise) averaged 
density p(r, t). The equations for the instantaneous microscopic density /)(r, t) are of little 
use because p{r,t) directly encodes the positions of all particles (modulo permutations); 
it is simpler in this case to follow the particle trajectories directly by integrating equation 
(HJ). But, depending on the situation one is trying to describe, it can be reasonable to 
consider either the coarse grained density p{r,t) or the ensemble averaged density p{r,t). 
In fact, their time evolution equations, (1211) and (fTSt respectively, are identical except for the 
additional noise term in ^TH . Both (approximate) equations are consistent with equilibrium 
thermodynamics. While the equilibrium value of p(r, t) (i.e. in the limit t oo), which we 
denote p(r), is correctly given by the stationary point of the Helmholtz free energy functional 
= n JHIIIQIIDj the time evolution of the coarse grained density p{r,t) also includes 
fluctuations about the equilibrium profile p{r). In other words p(r) is the most probable 
density distribution in iy('^^^[p]. For non-interacting particles, as pointed out in reference 
im, the deterministic DDFT in equation (fTSl) reduces to the (exact) drift-diffusion equation 
for Brownian particles in an external potential. These observations suggest that the DDFT 
provides a good approximation to the exact time evolution of the particle density, particularly 
when the system is either close to equilibrium or when the particles interact weakly. 

We believe that the origin of the debate as to whether the DDFT should be a stochastic 
or deterministic equation is a matter of confusion between the three types of densities defined 
in this paper. In particular, the distinction between the instantaneous microscopic density 
p(r, t) and the coarse grained density p{r, t) is of great importance, as noted also in reference 
ifTSl . The confusion is most apparent in two recent contributions to the field, references 013. 
In both papers, the Langevin and Fokker-Planck equation for p{r,t), namely Q and ©J 
respectively, are derived correctly, but then the microscopic instantaneous density p{r,t), 
which is a sum of delta functions, is confused with the coarse grained density p{r,t). In 
addition to this confusion there is, in references 013, a lack of a clear distinction between 
the functional 7i[/3] and the free energy functional T[p]. Bearing in mind the results from 
equilibrium density functional theory |ITJI3 this distinction should be evident immediately. 

The authors of reference |I3 argue that a DDFT has to be stochastic in nature. This 
misconception has its roots in reference where the difference between p(r, t) and p(r, t) 
is not recognized at all. This difference is further obscured by the fact that the time 
evolution equation for the spatially coarse grained density derived in |l6l| involves a free energy 
functional whose form is strikingly similar to 7i[/9] in Q. However, the crucial difference 
is that Q contains the bare pair interaction potential while the functional derived in |l6| 
contains the direct pair correlation function, instead of V{r), and therefore leads to the proper 
thermodynamic equilibrium state. We wish to emphasize that a DDFT for p(r, t), the density 
distribution averaged over all realizations of the stochastic noise, must be deterministic. 
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* Brownian dynamics for ■ 
{ri{t)}i=i...N equation 




Noise averaged 
1IIIIIIIIIIIIIIII 



Fokker-Planck equation for 
w{{ri}) equation (7) 



Evolution equation for 
p{r, t) equation (15) 



Fokker-Planck equation for 
l¥[p(r)] equation (9) 



i 



^Coarse grained 



Langevin equation for 
p(r, t) equation (21) 



Fokker-Planck equation for 
W[p{r)] equation (22) 



Figure 1. Flow chart outlining the various possible routes for deriving a DDFT for pairwise 
interacting Brownian particles. The starting point for all attempts to derive a DDFT is equation 
Q (thick dashed line box at the top). The end point is either a theory for the coarse grained 
density p(r, t) (Langevin or Fokker-Planck equation, thick solid outlined boxes at the bottom) 
or, for the ensemble averaged density p{r, t), a deterministic equation, (thick solid outlined 
box on right). The Langevin and Fokker-Planck equation for the microscopic density p{r, t) 
are only intermediate results (dotted boxes). Steps 1, 2, 5, 8, and 9 are exact (in circles) 
whereas steps 3, 4, 6, and 7 involve approximations (in square boxes). Ensemble averaged 
theories are located in the vertically striped area while coarse grained theories are located in 
the horizontally striped area. 



In summary, one can justify both a stochastic as well as a deterministic DDFT, depending 
on the quantities of interest. In this sense, there is a confusion rather than a controversy about 
the nature of DDFT. 

Acknowledgments 

We thank R. Evans for fruitful discussions on the present work and AJA acknowledges the 
support of EPSRC under grant number GR/S2863 1/01 . 

References 

[1] Evans R 1979 Adv. Phys. 28 143 

[2] Evans R 1992 in Fundamentals of Inhomogeneous Fluids, edited by D Henderson (Dekker: New York), 
chapter 3 

[3] Munakata T 1989 J. Phys. Soc. Japan 58 2434 

[4] Kirkpatrick T R and Thirumalai D / Phys. A: Math. Gen. 22 L149 



Time dependent density functional theory 



10 



[5] Dean D S 1996 /. Phys. A: Math. Gen. 29 L613 

[6] Kawasaki K 1994 Phyisica A 208 35 

[7] Frusawa H and Hayakawa R 2000 J. Phys. A: Math. Gen. 33 L155 

[8] Marconi U M B and Tarazona P 1999 J. Chem. Phys. 110 8032 

[9] Marconi U M B and Tarazona P 2000 J. Phys.: Condens. Matter 12 A413 

[10] Archer A J and Evans R 2004 submitted 

[11] Risken H 1984 The Fokker-Planck Equation, (Springer: Berlin) 

[12] Rauscher M, Mecke K and Griin G 2004 to be published soon 

[13] Zinn- Justin J 2002 Quantum Field Theory and Critical Phenomena, (Oxford University Press: Oxford), 
chapter 4 

[14] Goldenfeld N 1992 Lectures on Phase Transitions and the Renormalization Group, (Addison Wesley) 

[15] Kawasaki K 1998 J. Stat. Phys. 93 527 

[16] Dhont J K G 1996 / Chem. Phys. 105 51 12 

[17] Percus J K 1976 J. Stat. Phys. 15 505 

[18] Perez-Madrid A, Reguera D and Rubi J M 2002 / Phys.: Condens. Matter 14 1651 

[19] Reguera D and Reiss H 2004 J. Chem. Phys. 120 2558 



